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Abstract: Cosmic rays are hampered by the Moon and a deficit in its direction is expected (the so-called Moon shadow). 
The Moon shadow is an important tool to determine the performance of an air shower array. In fact, the displacement of 
the shadow center, due to the bending effect of the Geomagnetic field on the propagation of cosmic rays, allows to set 
the energy scale of the primary particles inducing the showers observed by the detector. The shape of the shadow permits 
to determine the detector point spread function. The position of the deficit at high energy allows evaluating its pointing 
accuracy. Here we present the observation of the cosmic ray Moon shadowing effect carried out by the ARGO-YBJ 
experiment (Yangbajing Cosmic Ray Laboratory, Tibet, RR. China, 4300 m a.s.l., 606 g/cm^ ) in the multi-TeV energy 
region with high statistical significance (70 standard deviations). By means of an accurate Monte Carlo simulation of the 
cosmic rays propagation in the Earth-Moon system we have studied the role of the Geomagnetic field and of the detector 
point spread function on the observed shadow. 
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1 Introduction 

The analysis of the Moon shadow observed by an air 
shower array may provide unique information on its perfor- 
mance. At high energies, the Moon shadow would be ob- 
served by an ideal detector as a 0.52° wide circular deficit 
of events, centered on the Moon position. The actual shape 
of the deficit as reconstructed by the detector allows the de- 
termination of the angular resolution while the position of 
the deficit allows the evaluation of the absolute pointing ac- 
curacy. In addition, charged particles are deflected by the 
Geomagnetic field (GMF) by an angle depending on the 
energy. As a consequence, the observation of the displace- 
ment of the Moon shadow at low rigidities can be used to 
determine the relation between the shower size and the pri- 
mary energy. The same shadowing effect can be observed 
in the direction of the Sun but the interpretation of this phe- 
nomenology is less straightforward [1]. Finally, the Moon 
shadow can be exploited to measure the antiproton content 
in the primary cosmic rays (CR). In fact, acting the Earth- 
Moon system as a magnetic spectrometer, paths of primary 
antiprotons are deflected in the opposite sense with respect 
to those of the protons in their way to the Earth. This effect 
has been used to set limits on the p flux at Te V energies not 
yet accessible to balloon/satellite experiments [2]. 

In this paper we present the observation of the cosmic ray 
Moon shadowing effect carried out by the ARGO-YBJ ex- 
periment during the period from July 2006 to November 
2010. We report on the angular resolution, the pointing 



accuracy and the energy scale calibration of the detector 
in the multi-TeV energy region. The results are compared 
with the predictions of a detailed simulation of cosmic ray 
propagation in the Earth-Moon system. 

2 The ARGO-YBJ experiment 

The detector is composed of a central carpet 74 x 78 m^, 
made of a single layer of Resistive Plate Chambers (RPCs) 
with about 93% of active area, enclosed by a guard ring 
partially (20%) insti-umented up to -- 1 00 x 11 m^ . The ap- 
paratus has a modular structure, the basic data acquisition 
element being a cluster (5.7x7.6 m-^), made of 12 RPCs 
(2.85 X 1.23 m^ each). Each chamber is read by 80 exter- 
nal strips of 6. 75x61. 80 cm-^ (the spatial pixels), logically 
organized in 1 independent pads of 55. 6x61. 8 cm^ which 
represent the time pixels of the detector. The read-out of 
18360 pads and 146880 strips are the experimental output 
of the detector [3]. The central carpet contains 130 clusters 
(hereafter, ARGO- 130) and the full detector is composed 
of 153 clusters for a total active surface of ^6700 m^. 

The whole system, in smooth data taking since July 2006 
firstly with ARGO- 130, is in stable data taking with the 
full apparatus of 153 clusters since November 2007 with 
the trigger condition Nt^jg = 20 and a duty cycle >85%. 
The trigger rate is ^-^3.5 kHz with a dead time of 4%. 

The analysis reported in this paper refers to events collected 
after the following selections: (1) more than 20 strips on 
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the ARGO-130 carpet; (2) zenith angle of the shower ar- 
rival direction less than 50°; (3) reconstructed core posi- 
tion inside a 150x150 m^ area centered on the detector 
According to the simulation, the median energy of the se- 
lected protons is E50 «1.8 TeV (mode energy «0.7 TeV). 



3 Monte Carlo simulation 

A detailed Monte Carlo (MC) simulation has been per- 
formed in order to propagate the CRs in the Earth-Moon 
system [4]. The simulation is based on the real data acqui- 
sition time. The Moon position has been computed at fixed 
times, starting from July 2006 up to November 2010. Such 
instants are 30 seconds apart from each other. For each 
time, after checking that the data acquisition was really run- 
ning and the Moon was in the field of view, primaries were 
generated with arrival directions sampled within the Moon 
disc according to the effective exposure time. 

After accounting for the arrival direction correction due to 
the magnetic bending effect, the air shower development 
in the atmosphere has been generated by the CORSIKA 
V6.500 code with QGSJET/GHEISHA models [5]. CR 
spectra have been simulated in the energy range from 10 
GeV to 1 PeV following the compilation results given in 
[6, 7]. About 10* showers have been sampled in the zenith 
angle interval 0-60 degrees. The experimental conditions 
(trigger logic, time resolution, electronic noises, etc.) have 
been reproduced by a GEANT4-based code [8]. 

The trajectory of a primary cosmic ray (energy E, charge 
Z) crossing the GMF is observed by a detector placed 
at Yangbajing bent along the East-West (EW) direction, 
whereas no deviation is expected along the North-South 
(NS) one (see Fig. 4, upper plot). To perform an evaluation 
of the bending effect, we adopted the Tsyganenko-IGRF 
model [9], which takes into account both internal and ex- 
ternal magnetospheric sources by using data available from 
spacecraft missions. To a first approximation, the amount 

of the EW shift can be written as: Aa ~ -1.58° jt \u - 

E[TeV] 

The sign is set according to the usual way to represent the 
EW projection of the Moon maps. This equation can be 
easily derived by assuming that the GMF is due to a pure 
static dipole lying in the centre of the Earth and is valid 
for nearly vertical primaries with energy greater than a few 
TeV. 

The contribution of different cosmic ray primaries (pro- 
tons. Helium and CNO group) to the Moon shadow deficit 
is shown in Fig. 1. Events contained in an angular band 
parallel to the EW axis and centered on the observed Moon 
position, compatible with the multiplicity-dependent angu- 
lar resolution, are used. 

The angular width of the Moon contributes to the spread of 
the signal, therefore we must disentangle this effect in mea- 
suring the detector angular resolution. Assuming a Gaus- 
sian PSF with variance a\, the width of the observed signal 
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Figure 1: Simulated deficit counts around the Moon pro- 
jected to the EW axis for N strip >100. The contribution 
of different primaries to the Moon shadow deficit can be 
appreciated. The CNO component has been multiplied by 
a factors. 



results: 



RMS = ae\ I 



0.52° 
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The contribution of the Moon size to the RMS is dominant 
when ae is low, i.e. at high particle multiplicities. For 
instance, the difference between RMS and <jg is 20% if 
ae = 0.2°, less than 5% if ae >0.4°, and only 1.7% if ae 
= 0.7°. The effect of the detector angular resolution along 
the EW projection of the Moon shadow deficit determines 
not only the smearing, but also a further displacement of 
the signal peak due to the folding with the asymmetrical 
deflection induced by the GMF. The West tail of the shifted 
signal, indeed, has a larger weight compared with the sharp 
East edge and tends to pull the signal in its direction. 

4 Results 

In Fig. 2 the significance map of the Moon region observed 
with data recorded in the period July 2006 - November 
2010 is shown for events with fired strips N strip > 100. 
The opening angle ip used in the smoothing procedure is 
1°. The statistical significance of the maximum deficit is 
about 70 standard deviations. The ARGO-YBJ experiment 
is observing the Moon shadow with a significance of about 
9 standard deviations per month. 



4.1 Angular resolution 

The PSF of the detector, studied in the NS projection not 
affected by the GMF, is Gaussian for N strip >200, while 
for lower multiplicities is better described for both MC 
and data with an additional Gaussian, which contributes 
for about 20%. For these events the angular resolution is 
calculated as the weighted sum of the ag of each Gaus- 
sian. In Fig. 3 the angular resolution measured along the 
NS direction is compared to MC predictions as a function 
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Figure 2: Significance map of the Moon region for events 
with N strip >100, observed by the ARGO-YBJ experi- 
ment in the period July 2006 - November 2010. The coor- 
dinates are R.A. and DEC. centered on the Moon position. 
The color scale gives the statistical significance in terms of 
standard deviations. 
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Figure 3: Measured angular resolution of the ARGO-YBJ 
detector (upward black triangles) compared to expectations 
from MC simulation (downward red triangles) as a function 
of the particle multiplicity. The multiplicity bins are shown 
by the horizontal bars. 



of the particle multiplicity, i.e. the number of fired strips 
^ strip on ARGO-130. The values are in fair agreement 
showing that the ARGO-YBJ experiment is able to recon- 
struct events starting from only 20 particles spread on an 
area ^6000 m^ large with an angular resolution better than 
1.6°. 

4.2 Absolute energy scale calibration 

In order to calibrate the absolute energy scale of CRs ob- 
served by the ARGO-YBJ detector we can use the GMF 
as a magnetic spectrometer In fact, the westward dis- 
placement of CRs by an angle inversely proportional to 
their energy provides a direct check of the relation be- 
tween the shower size and the primary energy. In Fig. 4 
the displacements of the Moon shadow in both NS (up- 
per plot) and EW (lower plot) directions as a function of 



1.25 1.76 2.30 3.55 6.05 12.2 26.0 E^JJeVyZ 




multiplicity 



Figure 4: Measured displacements of the Moon shadow 
(upper plot: NS, lower plot: EW) as a function of multi- 
plicity (black squares). The data are compared to MC ex- 
pectations (red circles). In the upper plot the solid line is 
fitted to the data. In the lower plot the solid curve is fitted 
to the MC events and the dashed curves show the ±10% 
deviation from the solid one. The energy scale refers to the 
rigidity (TeV/Z) associated to the median energy in each 
multiplicity bin (shown by the horizontal bars). 

the particle multiplicity are shown. The upper scale refers 
to the rigidity (TeV/Z) associated to the median energy 
in each multiplicity bin. The observed shift is compared 
to the results of the MC simulation of CR propagation in 
the Earth-Moon system. A shift of (0.19±0.02)° towards 
North can be observed. This displacement is independent 
of the multiplicity. Concerning the EW direction, the good 
agreement between data and simulation allows the attri- 
bution of this displacement to the combined effect of the 
detector PSF and the GMF. Therefore, the energy scale 
can be fixed in the multiplicity range 20-2000 particles, 
where the Moon shadow is moving under the bending ef- 
fect of the GMF. The MC results are fitted by the func- 
tion Aa = K{Nstrip)^, with k = -10.17 and A = -0.63, 
shown by the solid curve in Fig. 4, lower plot. To estimate 
the possible shift in particle multiplicity between data and 
simulation, as shown by the dotted curves in Fig. 4, the 
experimental data are fitted by the same function but with a 
multiplicity shift term: -10.17[(1 - ARn)Nstrip]~°-^^, as 
described in [10]. The parameter Ai?„ is the multiplic- 
ity shift ratio, resulting in Ai?„ = (+4 ± 7)%. Finally, 
the conversion from Ai?„ to the energy shift ratio ARe 
is performed. To determine the relationship between Ai?„ 
and ABe, and to check that this method is sensitive to en- 
ergy, six MC event samples in which the energy of the pri- 
mary particles is systematically shifted event by event in 
the Moon shadow simulation are calculated [10]. These six 
ARe samples correspond to ±20%, ±15% and ±8%. Fi- 
nally, by assuming a linear dependence, the relation Ai?„= 
(-0.91 ± Q.16)-ARe is obtained. Hence, the systematic 
uncertainty in the absolute energy scale ARe is estimated 
to be (+5 ± 8)%, where the error is the statistical one. 

Two systematic uncertainties may affect this analysis, the 
first related to the assumed primary CR composition and 



AUTHOR et al. PAPER SHORT TITLE 



the second to the use of different hadronic interaction mod- 
els. We evaluated both contributions following [10]: <Jchem 

= 7%, <Jhadr = 12%. 

Finally, the absolute energy scale uncertainty in the 
ARGO-YBJ experiment is estimated to be smaller than 

13% in the en- 
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ergy range from 1 to 30 (TeV/Z), where the Moon shadow 
is shifted from its position due to the effect of the GMF. 

4.3 Long-term stability of the detector 
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Figure 5: Upper panel: Displacement of the Moon shadow 
from the apparent center in the NS direction as a function 
of the observation time. Middle panel: Displacement of the 
Moon shadow from the apparent center in the EW direction 
as a function of the observation time. Lower panel: The 
ratio of the observed deficit count to the expected one as a 
function of the observation time. The plots refer to events 
with a multiplicity '^ strip >100. 

The stability of the detector performance as far as the point- 
ing accuracy and the angular resolution are concerned is 
crucial in 7-ray astronomy. Since November 2007 the full 
detector is in stable data taking with duty cycle >85%. 
Therefore, the stability of the ARGO-YBJ experiment has 
been checked by monitoring both the position of the Moon 
shadow, separately along R.A. and DEC. projections, and 
the amount of shadow deficit events in the period Novem- 
ber 2007 - November 2010, for each sidereal month and for 
events with '^ strip >100. 

The displacement of the center of the Moon shadow in the 
NS direction enables us to estimate the systematic error 
in pointing accuracy and its long-term stability aside from 
MC simulations, since the EW component of the GMF is 
almost zero at Yangbajing. The displacement of the shadow 
position from the Moon center in the NS direction is plotted 
in the upper panel of Fig. 5 as a function of the observa- 
tion time. Assuming a constant function, the best fit result 



(continuous line) shows that the Moon shadow is shifted 
towards North by (0.19±0.02)°. The RMS around this po- 
sition is 0.13°. 

In the middle plot the displacement along the EW direction 
is shown. The best fit result (continuous line) shows that 
the Moon shadow is shifted westward by (-0.36±0.02)°, in 
agreement with the MC expectations ((-0.35±0.07)°). The 
RMS around this Moon position is 0. 1 1°. 

The amount of CR deficit due to the Moon provides a good 
estimation of the size of the shadow, therefore of the angu- 
lar resolution. The expected number of deficit events "^def 

is approximatively given by Ndef — 1 — e^°'^'(^^ 
Nmoon, where R is the angular distance from the source, 
Nmoon is the number of events intercepted by the Moon and 
(J is the Gaussian width of the shadow. In the lower plot of 
Fig. 5 the ratio of the observed deficit count to the expected 
one is shown as a function of the observation time. The ex- 
pected deficit events have been calculated for (t=0.8°. The 
line shows the best fit result assuming a constant function: 
1.005±0.016. The RMS of the corresponding distribution 
is 0.11. 

5 Conclusions 

The Moon shadowing effect on cosmic rays has been ob- 
served by the ARGO-YBJ experiment in the multi-TeV en- 
ergy region with a statistical significance greater than 70 
standard deviations. 

We have estimated the primary energy of the detected 
showers by measuring the westward displacement as a 
function of the multiplicity, thus calibrating the relation be- 
tween shower size and CR energy. The systematic uncer- 
tainty in the absolute energy scale is evaluated to be less 
than 13% in the range from I to 30 (TeV/Z). 

The position of the Moon shadow measured with the 
ARGO-YBJ experiment turned out to be stable at a level of 
0.1° and the angular resolution stable at a level of 10%, on 
a monthly basis. These results make us confident about the 
detector stability in the long-term observation of gamma- 
ray sources. 
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